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Abstract 

The dynamic conductance of the heavy-fermion superconductor CeCoIns is measured by point- 
contact spectroscopy as a function of temperature from 60 K down to 400 mK. The contact between 
the gold tip and the single-crystal is shown to be in the Sharvin limit with the enhanced sub-gap 
conductance arising from Andreev reflection. The zero-bias conductance data are best fit using the 
extended Blonder-Tinkham-Klapwijk model with a d-wave superconducting order parameter. A fit 
to the full conductance curve at 400 mK indicates strong coupling {2A(0)/kBTc = 4.6) and quanti- 
fies the suppressed Andreev reflection signal, which is a signature of normal- metal/heavy- fermion 
superconductor junctions. Possibilites for theoretical modeling to account for the suppressed An- 
dreev conductance are suggested. 

PACS numbers: 74.50. +r, 74.45. +c, 74.70.Tx, 74.20.Rp 
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Among the heavy-fermion superconductors (HFSs), the recently discovered homologous 
layered family of CeMIns (M = Co, Rh, Ir) have been drawing particular attention because 
of the richness of the underlying physics (Ref. and references therein), including the 
coexistence of superconductivity and antiferromagnetism and the possible existence of a 
quantum critical point. The compound CeCoIus has been investigated most extensively 
because the TJ_2.3K) is attainable, and several experimental observations indicate 

either dx2_y2- or d^y- 4] wave pairing symmetry (PS). GoU et al. |5j have reported 

point-contact spectroscopy (PCS) data but without clear spectroscopic evidence of the PS. 

PCS has been widely adopted to investigate the PS of superconductors including HFSs 
0, ~\. Blonder, Tinkham, and Klapwijk (BTK) formulated a theoretical model for elec- 
tronic transport across a normal- metal /conventional superconductor (N/S) interface with 
an arbitrary barrier strength (Zefr) ^8]. A quasi-particle (QP) with energy lower than the 
superconducting gap energy (A), injected from the normal side, cannot enter the super- 
conductor as a single particle, but is retro-refiected as a quasi-hole, in order to conserve 
energy, momentum, and charge. This quantum mechanical phenomenon, called Andreev 

n 

reflection (AR) [9], is effectively a scattering from the superconducting pair potential. In an 
N/S contact with Z^^ = 0, the zero-bias conductance (ZBC) is predicted and observed to 
be twice the normal-state value 8]. Any mismatch in the Fermi surface parameters acts as 
an effective barrier, thus, reducing the probability of AR. In semiconductor-superconductor 
interfaces, results reported for Si-, GaAs-, InGaAs-, and InAs-based junctions with Nb 
counter-electrodes jlO] could be accounted for using the standard formula given by Blonder 
et al. Zes = [(1 — rY I'^r + ZlY^'^, where r = vfn/vfs, the ratio of the Fermi velocities 
of the two electrodes (note Zf,s remains invariant for r ^/r), and Zq represents the con- 
tribution of a dielectric barrier. However, since it is not possible to separate the effects of 
an impurity- or disorder-induced barrier (Zq) at the interface from that of the Fermi surface 
mismatch in these systems, the accuracy of Z^s remains inconclusive. In the case of an 
N/HFS interface, the BTK theory predicts that the conductance curve lies in the extreme 
tunneling regime {Z^s > 5) because of the large mismatch in Fermi velocities. However, 
it is common to observe AR-like enhancement of the sub-gap conductance (ESGC), albeit 
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suppressed in magnitude [ll|, [12, ll3|, [IJ, ll5|, ll6| . Deutscher and Nozieres [17[ address this 



inconsistency by proposing that the Fermi velocities entering in the ratio r are without the 
mass enhancement factor. 
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FIG. 1: Dynamic conductance spectra of a Au/CeCoIns point contact between 60 K and 400 niK. 
Curves are shifted vertically by 0.05 for clarity. 

In this Letter, we report PCS data with the tip parallel to the c-axis of the single crystal 
CeCoIns as a function of temperature from 60 K down to 400 mK. Electrochemically etched 
Au tips are used as counter electrodes. The tip-sample distance is adjusted electro mechan- 
ically |3l and the dynamic conductance (dl/dV) spectra are obtained by the standard 
four-probe lock-in technique. 

Figure 1 shows dl/dV versus V data normalized by the conductance at -2 mV. An 
asymmetry is seen to develop, starting between 40 K and 50 K. This is reminiscent of the 
emergence of a coherent heavy-fermion liquid at ~ 45 K |^19]. It becomes more pronounced 
with decreasing temperature down to 2.6 K, below which the background conductance re- 
mains almost the same. This behavior is attributable to the saturated relative weight of the 
coherent phase below ~ 2 K To facilitate analyses, we factor out the asymmetric part 
of the conductance data in Fig. 1 using the data taken at 2.6 K, obtaining fully symmetrized 
curves, as shown in Fig. 2. Two notable features in the 400 mK curve are the nearly fiat 
region near zero-bias and the ZBC enhancement (due to AR) of only about 13 %, suppressed 
heavily 13, 3| compared to those of N/S contacts with small Z^^. 

In order to check the nature of the contact, we estimate the contact size using Wexler's 
formula |2(i] and the high-bias {eV ^ A) resistance {Rn ~ 1.1^2), obtaining an upper 
hmit of 460 A. The elastic electron mean free path (EMFP), l^u is evaluated to be 810 
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A at Tc from the thermal conductivity (k) measurements for CeCoIns 21 1. We extend 



the calculation down to lower temperatures using the same data j21| and thermodynamic 



relations for low-energy QPs in a d-wave superconductor |22|. Namely, k/T oc p„r, pn oc 
T, T = lei/vp, .'. lei c<. k/T"^, whcrc T is the temperature, p„ the normal QP density, and 
r the QP lifetime. This results in l^i increasing exponentially with decreasing temperature, 
ranging 4-5 fim at 400 mK, nearly two orders of magnitude larger than the contact size. We 



also estimate the inelastic EMFP based on the microwave conductivity data 2^, obtaining 
a lower limit of 6500 A at 400 mK. Therefore, we conclude that the measured contact is 
truly in the ballistic (or Sharvin) limit at low temperatures, even if we take into account 
the possibility of reduced EMFPs in a point contact compared to those in a bulk. This 
extreme cleanness of CeCoIns together with Pauli-limited upper critical field enables the 
long-standing Fulde-Ferrell-Larkin-Ovchinikov phase transition to be observed for the first 
time in this material. Consequently, the arguments proposed by Gloos et al. 2^, 
attributing the suppressed ESGC to the non-ballistic nature of the contact, are not valid for 
our measurements. 

Detailed information of the electronic structure may be obtained by fitting the conduc- 
;ance curves using the extended BTK (EBTK) model formulated by Tanaka and Kashiwaya 
3|. Since the QP current is injected along the c-axis of CeCoIus, we integrate the kernel 
over the full half of the momentum space: 



ljl''d<P§dej^^dE{f{E-eV)-f{E)}as{E,4>)coseM 
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where as{E,(j)) = |i+^2(i^a2)|2 % A = — , (Tn = j:^, and Z = For d-wave 

PS, A(T, 0) = A(T)cos2(/). In order to incorporate a QP lifetime broadening, we replace 
E = E' — iP, where T = h/r is the QP scattering rate, and take the real part of the kernel 

El- 
Fitting is performed by numerical integration of Eq. ((T)) using Z^^A, and P as parameters. 

Although there are strong evidences for rf-wave PS in CeCoIus yiyil^y], s-wave fitting is 

also performed for completeness. Owing to the constancy of Rn, it is reasonable to set Zes 

to a constant 0.346, because the measured data are AR-like, thereby requiring small value 

of Zes 1^ ■ We also find that varying Z^s does not result in any better fits for either the s- or 

the (i-wave model. The best fit curves using the s-wave model are displayed as solid lines in 
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Fig. 2(a). The central flat region in low temperature curves can be reproduced by adjusting 
A and F. Note the calculated curve flts the data well near Tc, but deviates around the 
gap edge with decreasing temperature. Since the contact is ballistic, the local Joule heating 
effect is ruled out 0, as an origin of this deviation. The fltting parameters, A and F, 
are plotted in Fig. 2(b) as a function of temperature. From the value of A extrapolated 
to zero temperature, A(0) = 404 fieV, we obtain the ratio 2A{0)/kBTc = 4.08, which 
indicates strong coupling in accord with other experiments l29|l. We note F decreases with 

nn 

increasing temperature which is unphysical and in contrast to usual observations |27l . I30l | 
that F increases with increasing temperature. We attribute this behavior to the failure of the 
s-wave BTK model. In the case of (i-wave model, because of the contact conflguration, the 
PSs d^2_y2 and d^y are not distinguishable. The rf-wave flt curve for 400 mK is displayed as a 
solid line in Fig. 3, together with the s-wave flt curve and the measured data. We point out 
that the shallow dip seen around -1.2 mV in the data is not an intrinsic feature indicative of 
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FIG. 2: (a) Fitted curves (solid lines) to normalized conductance spectra (symbols) using s-wave 
BTK model. Curves are shifted for clarity, (b) Fitting parameters, A and F, along with BCS 
energy gap (solid line) and F (dashed line) used for ZBC fitting in Fig. 4. 
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FIG. 3: Comparison of best fit curves to the conductance at 400 mK. Dotted(solid) curve is for 
s(d)-wave EBTK fitting. 

the local heating effect 0, [2^, but an artifact caused in the normalization process due to an 
imperfect match of the background conductance. The fitting parameters are, Z^s = 0.365, 
r = 218 ^eV, and A = 460 fieV, which gives the ratio 2A/k-BTc = 4.64, implying again 
strong coupling 29]. The ci-wave model gives a better fit than the s-wave model, showing 
less deviation (albeit still substantial) above the gap edge and reproducing a slight dip-peak 
feature near the zero-bias. However, we find that the rf-wave model fits to the temperature 
dependence data only with decreasing F with increasing temperature, which is unphysical, 
similarly to the s-wave fit in Fig. 2(b). We interpret this as a failure of the EBTK model 
to explain the reduction in both the energy and the ESGC in our conductance curves. 



The conductance at the zero-bias would be least affected by any local heating effect 
The ZBC vs. temperature curve is fit to the s-wave model using A(0) = 349 fieV and Zcs = 
0.346, as shown in Fig. 4. The best fit is obtained only with decreasing F with increasing 
temperature, as plotted in Fig. 2(b): F(t) = 0.86A(0)(1 - tV3), where t = T/T^. This 
unphysical fitting parameter implies again a breakdown of the s-wave BTK model. For the d- 



wave model, the fitting parameters are Zcs = 0.365, A(T) = 2.35^6^0 tanh(2.06-y/Tc/T — 1), 
and F = 218 fieV. This constant F in the d-wave fit is not unphysical, in contrast to the 
s-wave. Thus, we argue that d-wave is a more likely PS, consistent with the literature 

mm 

However, the c?-wave EBTK model does not fully account for our data, as seen above. In 
the following, we consider possible origins of this failure. First, we note that our AR-like 



conductance spectra, similar to the PCS data on other N/HFS point contacts 
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IJ, ll^ll^, are in agreement with Deutscher and Nozieres' arguments [12|, in the sense that 
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FIG. 4: Comparison of best fit curves to the ZBC vs. temperature data. Dotted(solid) curve is 
for s(d)-wave fitting. 

they all fall to the weak-barrier region. A common observation in these PCS measurements 
is that the ESGC, which is 100% in an ideal N/S contact, is suppressed heavily by an order 
of magnitude. To our best knowledge, our PC spectra are one of the cleanest examples for 
the suppressed AR conductance (13%) with a complete data set over a wide temperature 
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are ruled out, since the ESGC 



range. In our case, proximity j31|] and pressure effects 

is observed to accompany the superconductingtransition. Moreover, the ballistic nature of 
;he contact excludes the local heating effect jo, and the dominant Maxwell resistance 
2^. Therefore, we claim that there must be intrinsic origins causing AR conductance in 

N/HFS contacts to be reduced severely. 



Golubov and Tafuri j33| consider a breakdown of the Andreev approximation (retro- 
reflectivity) when A/Ep {Ep is the Fermi energy) is not negligible and/or the electrodes 
have layered structures like the high-Tc cuprate superconductors. Since HFSs are known to 



so_be affected by this non-retrorefiectivity, 
consider mismatches in Fermi velocities 



have relatively small Fermi energies, they may 
resulting into a reduced ESGC. Mortensen et al. 
and momenta, showing reduced ESGC. However, it is the SGC normalized by a high-bias 
conductance [eV ^ A) that is suppressed, not the one normalized by the normal state 
conductance. Furthermore, considering the fact that the full conductance curve seems to 
be reduced in both the energy and the conductance, it is unlikely that the measured data 
can be accounted for with only mismatches in Fermi surface parameters In the 

case of inhomogeneous superconductors as in the two-fluid model jl^, it may be justifiable 
to set different barrier strengths for electrons and holes 3^ in each phase. In this case. 
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ihe asymmetry in the un- normalized data can also be accounted for. Anders and Gloos 
35^ put forward a theory, which accounts for both the reduced ESGC and the energy gap. 
They incorporate a strongly energy dependent F in HFSs to explain the strongly suppressed 
AR signal and attribute the reduced gap to the renormalization effect due to the strongly 
reduced QP spectral weight. A more rigorous and detailed theoretical modeling is needed 
to investigate these possibilities. In addition, the relevant time scales for a two-particle 
AR process, in contrast with a single particle tunneling process [2^, should be taken into 
consideration as well as the directionality of charge transport due to quasi-two dimensional 
nature of the Fermi surfaces and the effect of non- Fermi liquid nature of CeCoIns j^^l- 
In conclusion, PCS has been performed on the HFS CeCoIus between 60 K and 400 mK. 
Fits using EBTK models to the conductance spectra indicate that the superconducting PS 
is more likely rf-wave, in agreement with the literature. The heavily suppressed ESGC, a 
common observation in N/HFS point contacts, is not explained by existing models. Various 
possibihtes are suggested for future theoretical works. 
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